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Abstract 

It is common in Physiological Optics to fit the 
corneal and the lens surfaces to conic-based 
surfaces (usually ellipse-based surfaces), ob- 
taining their characteristic radius of curvature 
and asphericity. Here we show that the vari- 
ation in radius and asphericity due to experi- 
mental noise is strongly correlated. This corre- 
lation is seen both in experimental data of the 
corneal topographer Pentacam and in simula- 
tions. We also show that the effect is a char- 
acteristic of the geometry of ellipses, and not 
restricted to any experimental device or fitting 
procedure. 

1 Introduction 

Conic curves are the resultant of the intersec- 
tion of a cone and a plane, and have simple 
mathematical expressions [Ij. Ellipses are the 
most relevant conies for Physiological Optics, 
because both corneal and lens surfaces have 
usually elliptical profile. The general equation 
of an ellipse is 
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where (xq, yo) is the center of the ellipse, and a, 
b are the two semiaxes, Figure [TJ Frequently 
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(for example when we use the ellipse to de- 
scribe the profile of a lens) we are not inter- 
ested on the complete ellipse, but only in a part 
of it, usually near one apex (for example, thick 
region in Figured]). It is convenient to describe 
the apical region of an ellipse in terms of two 
parameters, radius of curvature and aspheric- 
ity [1]. The radius of curvature is the radius 
of the circumference that best approximate the 
ellipse in a very small region around the apex. 
The asphericity is a parameter that measures 
how the ellipse deviates from that circumfer- 
ence. We can define these two parameters in 
terms of the two semiaxes of the ellipse, as 



a 



(2) 
(3) 



where R is the radius of curvature and Q is the 
asphericit}Q, and both of them are referred to 
the superior apex of the ellipse (thick region in 
Figure [H) . 

Many surfaces, such as lens surfaces, are con- 
veniently described in terms of conic-based sur- 
faces (most commonly, ellipse-based surfaces). 
One example of these kind of surfaces is the 
ellipsoid (Figure [2]) , whose general equation is 
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* There are several different parameters for the de- 
scription of asphericity,all of them essentially equivalent 
[T], although see [2]. 
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Figure 1: Ellipse centered at {xo,yo) with semiaxes a, 
b (described by Equation [l}. 

However, this is not the only way of building a 
surface based on ellipses, and frequently used 
alternatives are, for example, biconical surfaces 
[31 HIE] and conicoids [H]. Again, frequently we 
are interested on the region around one apex of 
the surface. The geometry of this region is best 
described in terms of the radius and asphericity 
of the representative meridians. 

Ellipse-based surfaces parametrized in terms 
of radius of curvature and asphericity are 
widely used, for example to describe the sur- 
face of the cornea [D El H [3 El [71 [8] and the 
crystalline lens in the eye [9]. 

2 Fits to conic-based surfaces 

The problem that we will discuss here refers 
to the fitting of the experimental measurement 
of a surface to an conic-based surface, obtain- 
ing its representative radius of curvature and 
asphericity. All the fits described in this pa- 
per were performed by minimizing the squared 
distances in vertical direction (z axis, see Fig- 
ure[2]) between the experimental points and the 
fitted surface. 

This paper was motivated by observations 
on experimental measurements taken with the 
corneal topographer Pentacam (Oculus GmbH, 
Wetzlar, Germany). Therefore, we will focus 
on experimental datasets in the same format as 
those yielded by Pentacam. These data consist 




Figure 2: Computer-generated dataset. Points are on 
the surface of a perfect rotationally symmetric ellipsoid, 
with R — 6 mm and Q — 0.3. 

of a set of points in a 3-dimensional space, that 
are distributed on the X-Y plane on a regular 
square grid of side 100 /im, Figure [3 

3 Correlation between fitted radius 
of curvature and asphericity 

Figure [D shows a computer-generated dataset, 
following a perfect rotationally symmetric el- 
lipsoid, of radius Rreai = 6 mm and aspheric- 
ity Qreai = 0.3. If we fit these datapoints to 
an ellipsoid, we recover the right parameters, 
Rfit = 6 mm and Qfn = 0.3. Now, we add 
some noise to the datapoints, to mimic ex- 
perimental uncertainty. In particular, we add 
Gaussian noise with variance 0.01 mm to the 
z component of each point. If we take 1000 re- 
alizations of the noise, and fit the datapoints, 
we obtain some dispersion in the recovered pa- 
rameters, so that i?fit = 5.99 ± 0.02 mm and 
Qfit = 0.28±0.04 (mean it standard deviation). 
Figure [3l 

Now, we represent the fit's asphericity ver- 
sus the fit's radius of curvature. Figure [H We 
see a very strong and significant correlation, 
(r = 0.96, p < 10-l^ Q = -13.7 + 2.3R). We 
repeated this simulation with different ellip- 
soids and with different types of noise, get- 
ting the same correlation. Furthermore, we re- 
peated the simulation using fits to biconics [3] 
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Figure 3: Results of the fits to computer-generated 
datasets of rotationally symmetric ellipsoids with R = 
6 mm and Q = 0.3, with normal noise of variance 
0.01 mm added in the z direction, (a) Radius of curva- 
ture, (b) Asphericity. Both histograms are built from 
1000 simulations. 

instead of ellipsoids, getting the same result. In 
order to investigate whether this correlations 
were intrinsically motivated by the geometri- 
cal characteristics of ellipses, we performed the 
following calculation: We generated ellipse^ of 
radii of curvature ranging between 5900 and 
6200 /im, and asphericities ranging between 
0.1 and 0.5, and extending 3000 ^m around 
the apex. Then, we aligned them with re- 
spect to the ellipse with R = 6000 /xm and 
Q = 0.3, so that the squared deviations were 
minimal. And finally we computed the av- 
erage deviation between them and the ellipse 
with R = 6000 ^m and Q = 0.3. As was to 
be expected, this average deviation is higher 
the more different are the parameters with re- 
spect to the reference ones. The red contour in 
Figure H] encircles the region of the R-Q plane 
for which the average deviation is lower than 
0.5 nm. The good correspondence with the 
simulated data indicates that the correlation 
is due to this effect: The noise added to the 
surface causes the fit to find different param- 
eters, but with higher probability for ellipses 

^We perform this calculation with ellipses instead of 
ellipsoids (or any other ellipse-based surface) in order 
to distinguish the effects that are proper to elliptical 
geometry (and therefore are expected to reproduce in 
aU ellipse-based surfaces) from those that may be par- 
ticular of a certain type of ellipse-based surface. 
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Figure 4; Parameters of fits to ellipsoids of 1000 simu- 
lations with gaussian noise added to the surface. The 
red contour limits the region corresponding to ellipses 
whose average deviation with respect to the ellipse of 
R = 6000 fim and Q = 0.3 is less than 0.5 /im. 

that are most similar to the correct one. As 
these parameters for "similar ellipses" are lo- 
cated in a region of the R-Q plane that extends 
in a diagonal direction, the deviations due to 
the noise are correlated. 

We have observed the effect described here in 
measurements of the posterior corneal surface, 
taken with the topographer Pentacam (Ocu- 
lus) [8]. Figure [5^ shows the dispersion of re- 
peated measurements taken within few min- 
utes. Assuming that the eye does not change 
too much in such a short time, the observed 
dispersion is only attributable to experimental 
noise. Interestingly, the correlation also arises 
when these consecutive measurements are av- 
eraged, and we compute the change between 
measurements taken different days. Figure Eb- 
These measurements should be less influenced 
by measurement noise, due to the averaging. 
Thus, we might expect real changes in the 
cornea to be measured here, especially for eyes 
that correspond to patients that had under- 
gone refractive surgery (LASIK) between the 
two measurements. Although many of these 
changes are significant (see [8j ) , the correlation 
is still present. The correlation and the fact 
that overall dispersion is not much higher than 
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Figure 5: Results for experimental data of posterior 
corneal surface obtained with Pentacam. (a) Dis- 
persion of measurements taken consecutively, within 
few minutes. (b) Changes between measurements 
taken different days. Solid dots: Non-operated con- 
trols. Empty triangles: Patients, change one day after 
surgery. Empty circles: Patients, change one week af- 
ter surgery. Empty squares: Patients, one month after 
surgery. 

for the case of consecutive measurements indi- 
cates that part of the effect is due to measure- 
ment noise. However, note that some points 
deviate clearly more than for the case of con- 
secutive measurements. This may be indica- 
tive of some true change of the corneal ge- 
ometry, but the fact that the correlation still 
holds indicates that the change is probably 
very quite smallH 

4 Conclusions 

Dispersion in radius of curvature and aspheric- 
ity that are recovered by fitting an conic-based 
surface are strongly correlated. This is because 
the parameters that describe the most similar 
ellipses show the same correlation. This effect 

*We carried out a reanalysis of the data presented 
in [S], taking into account the correlation between R 
and Q. We find that, although the quantification of 
the changes may be biased (suggesting larger changes 
than those actually occuring), the distinction between 
significant and non-significant changes is essentially 
unchanged. Furthermore, as the correlation affects 
equally patients and control subjects (Figure [SJd) , the 
comparison between the two groups is valid. Thus, we 
find that the conclussins of [8] are essentially unmodi- 
fied by the reanalysis. 



should be taken into account when we try to 
describe geometrical of a surface, because the 
true change of geometry much be much smaller 
(or larger) than that apparent from the change 
in the two parameters. Also, this effect might 
be used to distinguish between real changes in 
the geometry of the surface and changes that 
are apparent from the fits, but are due to mea- 
surement noise. 
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